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The magnetic field anywhere inside a long solenoid with two symmetrically placed auxiliary 
windings is expressed in terms of zonal harmonics. Formulas are given that enable one to 
compute the size and position of the auxiliary windings that make the field most uniform 
at points in the neighborhood of the center of the solenoid. 

Inquiries from workers in the field of nuclear resonance have indicated their need of a 
region of uniform magnetic field of considerable volume. The field at any point in a solenoid 
can be calculated with an accuracy consistent with one's knowledge of its geometry and current 
but the uniformity of field inside ideal coils of reasonable length is not sufficient over the volume 
required to develop a satisfactory nuclear resonance signal at the 10 G level. 

The following analysis shows how the variable part of the field near the center of a solenoid 
can be nearly compensated by two symmetrically placed auxiliary windings. This method of 
compensation has been suggested and used by others 1 but a detailed analysis does not appear 
to have been published. 

The solenoid is treated as a current sheet made up of N turns of tape winding per unit 
length. The tape carries the current /, and the space between the turns of tape is assumed to be 
zero. The two equal and symmetrically placed auxiliary coils as shown in figure 1 with traces 
at T h 1\, and T[, T[ are treated as current filaments carrying the current I x each, where 
I X =N X I and N x is a disposable constant. The winding channels of the two auxiliary coils must 
each have linear dimensions of its axial section, say b and d, such that b 2 /a% and d 2 /a 2 t are so 
small compared with unity as to be negligible to the precision desired in the magnetic field. 
The radius a Y of each auxiliary coil and its axial position are to he determined so as to make 
the magnetic field most uniform near the origin. The equivalent filament that carries the 
current I x will have a trace that may be taken as the geometric center T x of the axial section 
of the winding channel. 

The points T x and r i\ are traces of a current filament of radius ai=r x sin 6 X in the plane 
x l =r 1 cos 0i. With spherical coordinates r, 6, and cylindrical coordinates x, p, </> the only com- 
ponent of the vector potential is A(r, d)=A (f) (r, 6) and this is independent of the longitude <j>. 
The spherical components of the magnetic field are 

Hr(r,e)= 7 4^0) e [A(r,e)sme] 

a(M)=72l P M(r,»)]. 

From these the cylindrical components are found by 

H x =H r cos 6 — lie sin 6 
H p =H r sin 6+ He cos 6. 



*Ketired. present address: 2520 Harrison Blvd., Ogden, Utah. 

i T. A. Meddle, Brit. J. Appl. Phys. 3, 95 (1952); B. Lewis, Brit. J. Appl. Phys. 1, 238 (1950); M. W. Garrett, J. Appl. Phys. 22 1091 (1951), 
M. R. Gardner, J. A. Jungerman. P. G. Liehtenstein. and C. G. Patten, Rev. Sri. Instr. 31, 929 (1960). 
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The vector potential at P(r, 0) due to a cgs unit circular current with trace at T x {r u 00 is 
given by 



A(r,e)=2irsm9sin 2 e l 51 (-) 

71 = 1 Vl/ 

=27rsin0sin 2 1 X(^\ 



r VPi (cos0)P;(cos r 



71(71+1) 



if 0<r<r 2 ; 



ri\ n+1 P f n (cos 6>)P; (cos 6>0 . 



n(n+l) 
where P„(u) is Legendre's polynomial and P' n (u) denotes -r~P n (u 



if r>r x , 



The effect of the two symmetrically placed equal filaments each with current I x is found 
from the first of these expansions, by superposition and noticing that in reckoning the effect of 
the filament on the left, the term P' n (cos 00 must be replaced by P' n (cos (?r— 00) =P' n ( — cos 00 
= (-l) n+1 P^(cos 00- 

The vector potential at P(r, 0) of both current filaments is therefore 



A il) (r, 0) = 4tt/ 1 sin sin 2 0i 2 ( - ) 

n=l Vi/ 



rY w_1 «.-i (cos 0)P^_ 1 (cos0O 



2n,(2w— 1) 



when 0<r<r! 



which applies therefore at points inside the solenoid near the origin. 
d 



By use of -j- [(1— u 2 )P' v (u)]= — v(v+\)P v (u) it is found that the field components are 



Hi l) =^I x 



sin 2 0i 



co / r \2: 



P 2n _, (cos 6)P' 2n _ l (cos 0,); 



H, (1) = -47ri lS in 2 ^ 



sin 6 - /r V" 1 PL-, (cos e)P^_, (cos t ) 






2rz— 1 



To obtain from these the expressions for the cylindrical components H x and H p it is neces- 
sary to make use of the two relations 



and 



We find 



{l-u 2 )P' v (u)+vuP v {u)^vP v ^{u) 

UP' V {U) -vP v {u) =P'„- 1 (U) . 



H?(r,e)- 



AtI x sin 2 0! 



[l+SGT 



-) P 2n (cose)PL+i (cos 00 



]i 



H™(r,ey- 



— 4ir/i sin 2 0! sin d » / r\ 2n PL (cos 0)PL-m (cos 00 



^i 



£(-) 



2n+l 



(1) 



The vector potential of the current / (cgs units) in the current sheet with N turns]per cm 
is found by integration to be 



co / r \2n p' 



1-1 co / r \\ 

A (s) (r, 6)=^NTr sin - cos a- sin 2 « 2 ( ") 



r\ 2n P L+i (cos d)PL (cos a ) 

2n(2n+l)(2n-\-2) 



] 



at a point P(r, 6) inside the current sheet. 
From this we find 



flj) ( r , 0) =47rA^/ Tcos a-sin 2 a 2 (0 



rV w P 2w (cos 0)P^ (cos ^""H 



2n 



] 



and 



^^^.aer *^^ 



(2) 
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For the case where the filament currents and solenoidal current circulate in the same 
sense around the x axis we place NJ=Ii in eq (1) and add eq (2). 

The resultant magnetic field inside the solenoid has the components 



tt / a \ a ATT f i n - 2 -2 ^ /r\ 2n P 2n (cos B)P 2n (cos a) 
H x (r, 6)=4:tNI < cos a +C sm 2 a-sin 2 a >; ( - 1 — — v 

rr / M A isjt • 2 • a ^ AY" K» (COS 6)P' 2n (COS a) A (7 \ 

H P (r,e)=4«NI sin* asms Jl(-) + ^ - (l-^J 



(-©} 



where 



C = 



Ni sin 2 Bi 
Nt x sin 2 a: 



62? 



•CD" 



PL (cos a) 



form=l,2,3, 



2nP' 2n +i (cos t ) 
Also we may place here for reference 

PJ(u)=3u 
«(«)= y (7« 2 -3) 



«(«i)=5 (5uJ-l) 
«(«i)=y (21ttJ-14«?+l) 



(3a) 
(3b) 

(4a) 
(4b) 



(5) 



For brevity, also let 



X = — > u=cos a, and w^cos 0i 



(6) 



where all three are positive (since X refers to the filament on the right in fig. 1). 

The variable part of H x is the series in (3a). If the first two terms corresponding to n — 1 
and n=2 are made to vanish, then the first two terms in the series (3b) for H p will also vanish. 
Therefore the approximately uniform resultant field near the origin will be in the axial direction. 

The two equations C 2 =C and C 4 = C are 



and 



2P' 3 (u 1 ) cN\\-u 2 Y 

yfiP'^u) _ N x (l-uXi 

4P' 5 (u 1 ) cN(l-u 2 )' 



At y o>0, 



(7) 
(8) 



Figure 1. Axial section of current sheet. 

c = the semi-diagonal of the section. 
a=the semi-aperture of its end-circle viewed from O. 
iV= number of turns per cm (tape-winding). 
7=the current circulating. 

The points (Ti, Ti (and) T\, T'\) are traces of the compensating 
filaments. 
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Dividing the members of eq (7) by the corresponding members of eq (8) gives 

X 2 



2 _2P' 2 (u)P' 6 (uQ 



or, in view of eq (5), 



r, / («;-l/3)*-4/ 63 

c \5(w 2 -3/7)(m?-1/5) w 

which must be real. We consider only the case of a relatively long solenoid in which 
u = cos a>Vr = 0.6547 or 0<a<49°6'. In order that X be real it is necessary that the value 
of u 1 ( = cos $i) to be found shall be such that 

(m?-1/3) 2 -4/63. 



(«?-l/5) 

The zeros of the numerator of this fraction are 



->o. 



|(l+-|)=0.5853 



V7 
and 

1 / 9. \ 

.08136. 



io-ih 



3V V7> 
Since in eq (9), U\ occurs in even powers only let 

Z=u\. (10) 

Equation (9) may be written 



n,_ Q / (g-0.5853)(g-0.08136) _ 

C ~ A "~^V 5(u 2 -f)(^-0.20) ' (11; 

If the value of z to be found makes X real, this will correspond to the case that was 
adopted above in placing I 1 = J r NiI so the current filaments circulate around the x axis in 
the same sense as the current / in the solenoid. If values of z are found that make X imagi- 
nary, a solution of the problem may still be found corresponding to I 1 = —Nil, which is 
equivalent to reversing the sign of C so that X will then be real. If X as given by eq (11) is 
inserted in eq (7), and the relations of eq (5) employed, then 

(l-s) W=f (l- S )( 3 -0.20)=^^^ 

3*cN(l-u 2 )u r (z-0.5853)(z-0.08136n 3 / 2 

2Nt ' L 5(u 2 -f) (g-0.20) J ' {U) 

Squaring eq (12) gives the following equation to determine z (which also applies to the 
case in which X is imaginary) : 

(l-g) 2 (g-0.20) 5 =5[(3-0.5853)(3-0.08136)] 3 (13) 

where 

Hence B is a positive constant if i£ 2 >f (or a < 49°6')- The numerical solution of eq (13) 
may be found without difficulty when B is given, that is, when the constants c, N, N 1} and 
u = cos a are given. 
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If we let 



F(z)>. 



_[(,3-0.5853)(2-0.08136)] 3/5 



and 



M=B~ 1/5 = 



(l-z) 2/5 

5( ?i 2_ 3)3/5 



rcNu(l^Y 



equation (13) is equivalent to 



M(z— 0.20) =F(z) wliere 2=cos 2 0,. 



(14) 



(15) 



(16) 



The roots of this equation between zero and one are the abscissas of the points of inter- 
section of the curve y=F(z) and the straight line y=M(z — 0.20). The curve F(z) could be 
plotted once for all, since the constants a and c of the solenoid enter only in the positive slope 
M of the straight line. The general nature of the curve y=F(z) is indicated in figure 2. A 
fairly accurate plot of F(z) and M(z — 0.20) versus z is helpful in getting a starting point for 
an accurate evaluation of the roots by successive approximation. If Z, is a root, the corre- 
sponding value of Ti is, by eq (12), 



KCKcidf^>-^-°- 2o > 



(17) 



which must be positive. There will be two roots of eq (16): Z, in the interval 0.")Sr ) 3<Z 1 <l, 
Z 2 in the interval 0.08136<Z 2 <0.20. The latter is inadmissible when Ii=+NJ. 



Figure 2. The curves y = F(z) and y = M(z-l/5) 
plotted as a function of z. 
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The root Zj gives the only value for cos 2 B x permitted when the currents in the filament 
circulate around the #-axis in the same sense as the solenoidal current /. If in the opposite 
sense, Z 2 would he indicated. When the solenoid and the value of the constant N x are arbitrarily 
chosen, this solution gives the best radius and position of the two current filaments in order to 
produce a uniform field in a sphere with center at 0. When the first two terms are th i is annulled 
in the series of eqs (3a) and (3b), the next term (n=3) will indicate the principal remaining 
nonuniformity of field near the center by comparison with the constant term (cos a-\-C sin 2 a) 
in eq (3a). 

The foregoing analysis has been applied to the compensation of a solenoid designated as 
NBS 58080. The constants of this solenoid at 25 °C are iV=9.99915 turns per cm, c = 51.9138 7 
cm, and cos a =0. 963064. For several assumed values of Ni the constants of the required 
auxiliary windings have been found and listed in table 1. For each value of N h the relative 
contribution of the compensating windings to the field at the center of the solenoid is shown 
in the column headed C sin 2 a/cos a. With the terms n = l and 2 removed from eqs (3a) and 
(3b), the value of r, at which the term n = 3 reaches one part per million (ppm) relative to the 
constant term (cos a+C sin 2 a) in (3a) will be found in the table under "compensated radius. " 
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Table 1. Constants of auxiliary windings for solenoid NBS 
58080 









Condition: I\ = 


+NJ 






















Co sin 2 a 


Com p. 


Ni 


M 


Z x 


X 


T\ 


n cos 0i 


ri sin 0i 




radius 


















cos a 












cm 


cm 


cm 




cm 


49 


0. 9946 


0. 7903 


0. 9423 


48.91 


43.5 


22 A 


0.022 


7.5 


16 


.6356 


.693 


.6938 


36.02 


30.0 


20.0 


.014 


5.8 


1 


. 2097 


.602 


.2805 


14.56 


11.3 


9.2 


.003 


3.1 


0.25 


.1204 


.592 


.1767 


9.17 


7.1 


5.9 


.001 


2.3 



Table 2. Constants of auxiliary windings for solenoid NBS 
58080 









Condition: Ii = 


= -XiI 






















Co sin 2 a 


Comp. 


Ni 


M 


z 2 


X 


n 


n cos 0i 


r\ sin 0i 




radius 


















COS a 












cm 


cm 


cm 




cm 


-49 


0. 9946 


0. 1139 


0. 8017 


41.6 


14. 05 


39.2 


-0. 108 


5.4 


-16 


.6356 


.1006 


.5823 


30.2 


9.59 


28.7 


-.059 


4.6 


-1 


.2097 


.0852 


. 2438 


12.7 


3.69 


12.1 


-. 0075 


2.7 


-0.25 


.1204 


.0829 


. 1547 


8.03 


2.31 


7.7 


-. 0029 


2.0 



Corresponding values for the condition I^—NJ have been found and are listed in table 2. 
There appears to be nothing in favor of the latter case, unless the relatively close spacing of 
the auxiliary windings is especially appropriate to an experiment. To give an idea of the 
expected improvement over the uncompensated solenoid, we calculate the term n=\ in the 
series of H x (eq (2)) relative to the constant term cos a, and find that it reaches 1 ppm at 
r=0.16 cm. 

We assume in what follows that the solenoid construction is accurate to 1 ppm. When 
deciding upon a design, one has to consider the accuracy required in the construction of the 
auxiliary coils in order to realize a given accuracy in the combination. If an accuracy of, say, 
1 ppm is desired and an experiment requires a uniform region whose volume has a radius of 
7.5 cm the required accuracy in the compensating coils (2Vi = 49) is about 1 part in 20,000, a 
rather severe requirement. On the other hand, a uniform region with 2 cm radius (iVi = 0.25) 
can be had for an accuracy of only 1 part in 1,000 in the construction of the auxiliary windings. 
Obviously it does not pay to compensate over a region larger than necessary. 

(Paper 69C1-185) 
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